Perturbation Spectra in the Warm Inflationary Scenario 
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We investigate the phenomenology of warm inflation and present generic results about the evolution 
of the inflaton and radiation fields. The general conditions required for warm inflation to take place 
are derived and discussed. A comprehensive approach is presented for the generation of thermally 
induced adiabatic and isocurvature perturbations and the amplitude of their spectra calculated. 
In addition we derive the ratio of tensor-to-scalar perturbations and find the spectral indices for 
adiabatic, isocurvature and tensor perturbations formed in the warm inflationary era. These results 
are presented in a simplified and compact approach that is generally applicable. Our results are 
illustrated by inflation models with polynomial and exponential potentials. We compare our analytic 
results against numerical models and find excellent agreement. Finally, presently available data 
is used to put constraints on warm inflation and we discuss how near-future observations may 
distinguish the warm inflationary scenario from standard supercooled inflation. The main observable 
difference is the different scalar-to-tensor ratio, and that the consistency relation between this and 
the slope of tensor perturbations does not hold for warm inflation. 
PACS number(s): 98.80.Cq, 98.80.-k, 98.80.Es 
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I. INTRODUCTION 

One of the central paradigms in cosmology is the idea 
of inflation, and its elegant solution to the problems of 
standard cosmology. In addition inflation provides an ex- 
planation for the origin of structure from quantum fluc- 
tuations frozen in by crossing the Hubble radius during 
the rapid expansion. However in the standard inflation- 
ary scenario (see e.g. ref the rapid expansion phase 
and subsequent reheating of the universe back to a radi- 
ation dominated model are treated as two separate pe- 
riods. The motivation for this was both simplicity and 
the belief that interactions with the field driving inflation 
would be minimal || . The switching on of field interac- 
tion dynamics at the end of inflation results in violent re- 
and pre-heating effects such as a resonant amplification 
which can lead to generic amplification of perturba- 
tions 

From the point of view of particle physics, interactions 
are inevitable, and they lead one to consider the effects 
on inflation of a phenomenological "friction" term de- 
scribing the decay of the inflaton field into other fields 
during inflation |lj,[ll|. One result of this is that den- 
sity perturbations are no longer generated via quantum 
fluctuations, but are dominated by the larger thermal 
fluctuations 

In this paper we consider the warm inflationary sce- 
nario |l5| , |l5| ] from the phenomenological point of view. 
Several types of phenomenological warm inflation mod- 
els exist in the literature p6|-p0[|. Warm inflation dy- 
namics also has been studied in quantum field theory M] 
and some of the phenomenological models can be derived 
from first principles [ ]l4| , p2| -p4| . Our purpose here is to 
determine generic observational features of the warm in- 
flation scenario, without referring to specific models. 



The paper is organised as follows. In Section II results 
are derive for evolution of the inflationary and radiation 
fields during the slow-roll phase of warm inflation. As- 
suming the dissipation rate dominates over the expan- 
sion rate, a relationship is found between the radiation 
energy density and the energy density of the inflation- 
ary field. The spectrum of adiabatic perturbations is 
calculated in Section IIIA and its spectral index in Sec- 
tion IIIB. Isocurvature perturbations are also generated 
during warm inflation and we present a new derivation 
of them in Section IIIC and of their spectral index in 
Section HID. Tensor perturbations are discussed in Sec- 
tion HIE. We compare our results with standard inflation 
in Section IIIF. The conditions for ending warm infla- 
tion are discussed in Section IV. Our results are illus- 
trated by the specific cases of inflation with polynomial 
and power-law potentials in Section V. Our analytic ex- 
pressions are tested against exact numerical models in 
Section VI, where we find good agreement between the 
two. In Section VII we discuss observational constraints 
on warm inflation and how it may be distinguished from 
the standard inflation scenario by future CMB experi- 
ments. Finally, our conclusions are presented in Section 
VIII. Two appendices are included in the paper. The 
first deals with a derivation of the thermal fluctuation of 
the inflationary field, while the second presents new re- 
sults for generic polynomial and exponential potentials. 
Throughout this paper, the conventions of |2q ] are fol- 
lowed. 



II. EVOLUTION OF THE INFLATIONARY AND 
RADIATION FIELDS 

The evolution of the inflaton field, </>, is given by 
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(3H + T)j> + V = 0, 



(1) 



where H = a/a is the Hubble parameter, a is the cos- 
mological expansion factor, and T is the dissipation co- 
efficient. V(cp) is the potential of the inflationary field. 
We assume a spatially flat universe. A dot denotes dif- 
ferentiation with time and a prime indicates d/dip. The 
energy-density of the radiation field, p 7 , obeys the trans- 
fer equation 



p\ + AHp^ = T<p 2 



(2) 



The equations are completed by the Friedmann equation 
for the evolution of the universe; 
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(V + P 7 ), 



(3) 



where m p is the Planck mass 



During an inflationary 
era the potential field dominates both the kinetic energy 
of the inflationary field and the energy density of the 
radiation, so the Friedmann equation can be reduced to 



(4) 



We assume the slow-roll condition, <j> 2 <C V(<fi), through- 
out the warm inflationary era, yielding 



V 



3H(l + r)' 



(5) 



where the ratio of the production rate of radiation to the 
expansion rate, r, is defined as 



r 



(6) 



and parameterises the importance of dissipation in the 
warm inflationary model. 

For the dissipation-dominated regime of warm infla- 
tion, the dissipation rate, T, is much greater then the 
expansion rate; r> 1. The evolution of the inflaton field 
during this phase is governed by 
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3rH' 



(7) 



Hence the slow-roll condition can be met even for rela- 
tively steep potentials, due to the damping effect of dis- 
sipation on the evolution of (f>. 

If we also assume that during warm inflation radiation 
production is quasi-stable (p' 7 <C r</> 2 ), equations (0) and 
& give 



lh = aT 4 = 3r ( - 
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(8) 



where a = g*it 2 /30 is the Stefan-Boltzmann constant and 
<?* is the number of degrees of freedom for the radiation 



field. In the standard cosmology « 100. Equation 
(H) shows that the energy density in the radiation field is 
controlled by the kinetic energy of the inflationary field as 
well as the dimensionless dissipation rate. As the kinetic 
energy of 4> is small during the slow-roll regime we require 
r to be large to generate substantial radiation energy- 
density. Equation (|§|) is not strictly true at the start 
and end of the warm inflationary phase, as the stable 
radiation production condition will be violated. 

In contrast the energy density of the inflaton during 
the slow-roll phase is dominated by its potential energy, 
giving 



P0 - V(<f>). 
The radiation density can be re-written as 
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(9) 



(10) 



during warm inflation, since in the slow-roll phase the 
kinetic energy of the inflaton is controlled by the slope of 
the inflationary potential. 

Introducing the dimensionless inflationary slow-roll pa- 
rameter 
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(11) 



we can relate the energy density of the radiation field to 
the energy-density of the inflaton by 



(12) 



We have verified by numerical modelling (Section VI) 
that equation ( |l2| ) holds very well during most of the 
warm inflationary phase. 

We assume that warm inflation begins from a radia- 
tion dominated era, and starts when the energy density 
in the inflationary field is equal to that in the radiation 
field. Initially the radiation energy is redshiftcd away 
by the quasi-exponential expansion but enters a quasi- 
static phase when radiation production stabilises during 
the dissipation-dominated phase. Warm inflation ends 
when the universe heats up to become radiation domi- 
nated again, at a time when 



P~< = P4>- 



(13) 



At this epoch the universe stops inflating and smoothly 
enters a radiation dominated big-bang phase. This is one 
of the most attractive features of warm inflation, and 
suggests that its predictions are robust in the sense that 
they are not altered during the smooth transition to the 
radiation-dominated epoch. The material components 
of the universe are created by the decay of either the re- 
maining inflationary field or the dominant radiation field. 

The approximate relation, equation (12), and equation 
( 13|) suggest that an approximate condition for warm in- 
flation to both begin and end is 
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£ = 2r, 



(14) 



III. THE PERTURBATION SPECTRA 



although we do not expect this relation to be accurate 
since the condition of stable radiation production is vio- 
lated at both epochs. In particular at the beginning of 
warm inflation the ambient radiation field is redshifted 
away during the sharp transition to warm inflation and 
the radiation density drops rapidly. Then at the end 
of the warm inflation phase, the production of radia- 
tion must end and we expect p 7 'SHp 7 w Tip 2 during 
the turnover regime going from inflation-domination to 
radiation-domination. 

In addition at both the beginning and end of warm 
inflation as the energy density of radiation is comparable 
to that of the inflationary potential, equation (||) should 
be used to describe the evolution of the universe rather 
than equation (^J), introducing another approximation. 

In Section VI we test the accuracy of our assumptions 
over the whole course of a warm inflationary phase by 
numerically integrating equations ([!]), (||) and (||) and 
find that for many models it appears that these approx- 
imations are accurate enough to be used to estimate the 
conditions over a large period of warm inflation, includ- 
ing the final few e-folds at the end of warm inflation. 

The condition for a warm inflationary phase can be 
summarised by 



£ < 2r. 



(15) 



It is straightforward to show that the same result can be 
derived from the requirement of accelerated expansion in 
the warm inflationary epoch. This requires that a > 
implying —H/H 2 < 1 or —4>H'{(j))/H 2 < 1. Using equa- 
tions (||) and (0) we again find the inequality of equation 
(|l5|). The slow-roll condition and the additional strong 
dissipation condition 



r > 1 



(16) 



then implies warm inflation. One immediate consequence 
is that the usual inflationary condition, e < 1, is no longer 
required as inflation (quasi-exponential expansion) can 
be maintained for steep potentials due to the damping of 
dissipation. 

Finally, one can differentiate equation (|7]) with respect 
to time and, following from the slow-roll condition tfi <C 
V , we find a constraint on the second derivative of the 
potential, 



where 



m' p V' 1 



(17) 



(18) 
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is the curvature slow-roll parameter for inflation. This 
should be regarded as an extra condition for sustained 
warm inflation, and is a generalisation of the usual in- 
flationary constraint, |?7| <C 1. Hence warm inflation can 
take place for a wider range of potentials than standard 
supercooled inflation. 



A. Adiabatic perturbations 

Density perturbations are calculated in the same way 
as for standard inflation Jl^,|l4| , via the equation 



r 2 (H 



(19) 



where 4> is given by equation (Q). This equation holds 
for a mixture of inflationary and radiation fields, in the 
quasi-static regime of radiation production. It also can be 
shown that this result holds for any scalar field interact- 
ing adiabatically with another field, in the quasi-static 
limit, with the equation of state p — up. In the pres- 
ence of non-adiabatic perturbations (see Section [II C| ) 
this relation does not necessarily hold. Specifically, if 
there non-adiabatic pressure perturbations, these can in- 
duce curvature perturbations to evolve on super-horizon 
scales [^6|. However, we expect these effects to be small 
for the scenario discussed here and equation ( |l9| ) can still 
be used for the order-of-magnitude estimates presented 
here. 

The amplitude of perturbations, Sjj, is related to the 
final dimensionless linear power spectrum of mass per- 
turbations by 



A 2 (fc) 
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T 2 (k)S 2 H (k) 



(20) 



where T{k) is the matter transfer function relating initial 
spectra to the final, post-recombination linear spectra. 

The fluctuations in the inflationary field are no longer 
generated by quantum fluctuations, but by thermal in- 
teractions with the radiation field; 



k F T 
27r r ' 



(21) 



where Uf = vTH = V 3rH > H is the freeze-out scale at 
which dissipation damps out the thermally excited fluc- 
tuations. The origin of the freeze-out wavenumber kp is 
explained in Appendix A based on the original deriva- 
tion in |lj|. Note that in the warm inflationary epoch, 
the freeze-out length scale is generally much smaller than 
the Hubble radius. As one final point, when T < H and 
T < H, the amplitude of thermal fluctuations is equal 
to the amplitude of quantum fluctuations generated by 
Hawking radiation. 

Combining equations (@), (@), © and (H) we nnd 
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(22) 



The temperature of the radiation field is related to the 
potential by equations (||), (||) and ([l2]) yielding 
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rpl 
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2ra 



(23) 



Expressing the amplitude of density perturbations in 
terms of inflationary parameters we find 
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(24) 



The ratio of thermal to quantum fluctuations is r 5 ^ 2 T/H. 
During strong dissipation, r» 1, while towards the end 
of the warm inflationary era T/H ~ rrip/V 1 ^ 4 1. 
Hence the thermal spectrum will dominate over the quan- 
tum spectrum of fluctuations during a warm inflationary 
phase. 



B. Spectral index of adiabatic perturbations 

The spectral index of the adiabatic scalar perturba- 
tions is given by 
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dink 



(25) 



where the logarithmic interval in wavenumber is related 
to the number of e-folds from the start of inflation, N(4>), 
and the scalar field by 
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d\nk(cf>) = -dN(4>) = j— < 

V 



(26) 



Hence the spectral index of adiabatic density perturba- 
tions in warm inflation is given by 



n s = l--(3e-2r 1 ). 



(27) 



As we require that r> lwe expect the deviation from 
scale invariance to be suppressed in the warm inflationary 
model, until e ~ r, at the end of the warm inflation phase. 
If e > (2/3)?7, or r\ is negative, the spectrum of adiabatic 
fluctuations will be red, with the spectrum suppressed on 
small scales, as in standard supercooled inflation which 
generally also produces slopes n s < 1. However, in the 
opposite case, when e < (2/3)?7, the spectrum of adia- 
batic fluctuations will be blue, in contrast with predic- 
tions from generic standard inflation models. With a blue 
spectrum, the interesting possibility exists of producing 
primordial black holes on small scales p7y. 



is undergoing thermal fluctuations relative to the number 
density of photons 



S 



x 



X 



(28) 



and we have assumed the X-field is the ambient radiation 
field excluding photons. 

Assuming that the isocurvature fluctuations are ther- 
mal in origin, their perturbation spectrum is 
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(29) 



This relation is true only for non-interacting fields, but 
again should hold for order-of-magnitude estimates. The 
ratio of isocurvature to adiabatic perturbations is 
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(30) 



so that when r ~3> 1 and r > e adiabatic perturbations 
will dominate. In practice the scalar modes will form a 
quadratic sum of adiabatic and isocurvature modes, 



S 2 S = % 
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(31) 



and which can be separated observationally since isocur- 
vature generated fluctuations of the CMB are out of 
phase with adiabatic perturbations. Another discrimi- 
nator is if super-horizon perturbations can be observed, 
for instance in the CMB [ p8[ . Since isocurvature modes 
do not contribute to curvature perturbations (in the co- 
moving gauge) the contribution from isocurvature modes 
should be absent from super-horizon perturbations, lead- 
ing to a drop in amplitude of the scalar perturbations at 
the horizon scale. Isocurvature perturbations have previ- 
ously been discussed in the context of warm inflation by 
p9[ , although no detailed calculation of their amplitude 
was given. 



D. Spectral index of isocurvature perturbations 

Following the derivation of the spectral index of adia- 
batic perturbation, the spectral index of the isocurvature 
perturbations, defined by 



dink 



(32) 



is 



C. Isocurvature perturbations 

Isocurvature perturbations also are generated in the 
warm inflationary era due to thermal fluctuations in the 
radiation field. These perturbations can be characterised 
by fluctuations in the entropy, Sx, where the species X 



4r 



{e + 2ri). 



(33) 



This differs from isocurvature fluctuations generated by 
quantum fluctuations, where n; so = — 2e. If r is large, 
the spectrum will again tend to be flat, until e ~ r. If 
e > — 2r\ the thermal isocurvature spectrum will be red 
as it is in standard supercooled inflation. 
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E. Tensor perturbations 



IV. CONDITIONS AT THE END OF INFLATION 



Tensor perturbations do not couple strongly to the 
thermal background and so gravitational waves are only 
generated by quantum fluctuations, as in standard su- 
percooled inflation. Their spectrum is given by 129] 



At = 



32V 



with spectral index 
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(34) 



(35) 



and the inequality — 4 < n g < arising from the warm 
inflation condition e < r. 

The ratio of tensor to scalar perturbations is 



4 _ ^ f2e 
9 ~ 6% ~ 3 \r* 
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(36) 



Hence, as we are in the strongly dissipative regime, where 
r 1 and s < r, we expect the contribution from grav- 
itational waves to be small unless the magnitude of the 
potential is large or r is order unity. 

The ratio of the amplitude of isocurvature to tensor 
perturbations scales as 
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(37) 



which will be greater than unity unless the scale of the 
potential is large. Hence we can expect that the inequal- 
ities 



Rg < Rh 



< 1 (38) 

in general will hold. We find that this is true for all the 
of the numerical models we consider in Sections V and 
VI. 



F. Comparison with Standard Supercooled Inflation 

A nice feature of warm inflation is that these results are 
significantly different from standard supercooled inflation 
case, where 
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n s = 1 — 6e + 2r), 
n„ = -2s, 



R a 



(39) 



Thus near future observations should be able to distin- 
guish between the two. An immediate consequence of our 
analysis is that there need be no consistency relationship 
between the scalar spectral index and the ratio of tensor 
to scalar perturbations, as in standard inflation. We dis- 
cuss the observational differences between the models in 
Section VII. 



At the end of the warm inflationary epoch, the dissipa- 
tion parameter, r, is related to the slow-roll parameter, e, 
by equation ([l4|). Substituting in equations (||), (||) and 
we form a differential equation which can be solved 
for V(4>) to find the magnitude of the potential field at 
the end of warm inflation, 



v{4> c 
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(40) 



up to a constant, but independent of the shape of the 
inflationary potential. This selects out the V ~ (p 4 po- 
tential as a special case where warm inflation never ends 
(see Appendix B). This is due to the energy density of 
the radiation field "tracking" the energy density in the 
inflationary field, p 7 oc p^, and so the condition for warm 
inflation to end is never met. In this case warm inflation 
continues until the slow-roll conditions are no longer met 
and the inflaton field decays away, which then ends infla- 
tion. For n > 4 and exponential potentials, the radiation 
density will always fall faster than the energy-density of 
the inflaton field (see Appendix B). In this case warm 
inflation will occur, in the sense that a sizable radiation 
component can persist for a considerable duration dur- 
ing inflation, well beyond the requirements for solving 
the problems of standard cosmology. However to end 
such a warm inflation regime, this behavior of the poten- 
tial must be local, and at some point the vacuum energy 
must be released at a sufficiently fast rate to allow the 
radiation to dominate. 



V. SPECIFIC CASES 

In order to understand how warm inflation differs from 
the standard supercooled inflationary scenario, it is use- 
ful to consider a few specific examples. Here we consider 
models of warm inflation with polynomial and exponen- 
tial potentials. 



A. Polynomial inflation 

We consider the simplest inflationary model 

V(<j>) = m 2 (/) 2 /2 (41) 

for a massive scalar field with mass m. With this poten- 
tial the slow-roll parameters are 



e = i] 
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The Hubble parameter is given by 
V 3 V m v 



(42) 



(43) 
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while the dissipation parameter scales like 



(44) 



From equation (|44| ) we see that the ratio T/m controls 
the evolution of the system. 

Solving equation (Q) we find that the inflaton field 
evolves as 



Compared with standard inflation, where JB| 



ct>{t) « 0o 



t, 



(45) 



(46) 



we see that during warm inflation the inflaton field decays 
due to dissipation into the radiation field. 

The energy densities of the inflationary potential is 
given by equation (|J) and ( f4f"| ) while the radiation field 
is given by 



' 3 / to to p 
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(47) 



The number of e-folds from the end of warm inflation is 
given by 

N(<f>) = -kik(<t>) 




fend 



(48) 



Figure 1 shows the evolution of the energy densities 
of the radiation and inflationary fields as a function of 
number of e-folds. The upper line is the inflation energy 
density and the lower line the radiation energy density. 
The dashed lines are the results of our analytic model, 
from equations (pill), ( f47| ) and (p0|), while the solid lines 



are the results of our numerical model (Section VI) of the 
evolution of the field through the warm inflationary era. 

From the condition that e — 2r at the end of warm 
inflation we find that the magnitude of the inflationary 
field at this time is 



(49) 



where the initial field must have |0o| > 0cnd- In stan- 
dard polynomial (n = 2) inflation this condition is 
1 01 end > Tip/V47r, so we find that for warm inflation 
this is weakened by the factor m/T. If T > to then 
warm inflation can begin well below the Planck mass. 





N=ln(a f /a) 

FIG. 1. The evolution of energy density of the inflationary 
and radiation fields as a function of number of e-folds from 
the end of warm inflation, N = ln(a//a), where a/ is the final 
expansion factor, for the inflationary potential V = m <j> /2. 
The energy-density of the inflationary field (upper curves) 
and the radiation field (lower curves) are both normalised 
to unity in units of the magnitude of the potential at the 
start of inflation, V(cj>o). The solid lines are the results of the 
analytic expressions equations (^?j) , while the dashed lines are 
the results of the numerical model discussed in Section VI 



The total number of e-folds from the end of warm in- 
flation is 




(50) 



As we require TV w 70 e-folds to comfortably solve the 
problems of standard cosmology this implies 



35 



(51) 



The requirement of strong dissipation also leads to a 
constraint on V. From equation (|44|) with equation ( |5l| ) 
we find that at the beginning of inflation 
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resulting in the requirement that T 3> 15m for warm 
inflation to commence in polynomial (n = 2) inflation. 

Calculating the spectrum of adiabatic perturbations 
yields 
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(53) 
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FIG. 2. The perturbation spectra of adiabatic scalar per- 
turbations (top line), isocurvature scalar perturbations (mid- 
dle line) and gravitational waves (bottom line) as a function 
of e-folds from the end of warm inflation, N, for the potential 
V — m 2 4> 2 /2. The model is normalised to have 8h = 2 x 10 -5 
at the 50 th e-fold. Solid lines are the results of our analytical 
model, while the dotted lines are from the numerical model. 



N=ln(a f /a) 

FIG. 3. The spectral indices of the adiabatic scalar pertur- 
bations, n s — 1, isocurvature scalar perturbations, ni so , and 
gravitational waves, n g , as a function of e-folds from the end 
of warm inflation, N, for the potential V = m 2 2 /2. The solid 
lines are for our analytic model, while dotted lines are from 
the numerical model. The spectral index of the isocurvature 
perturbations lies under that of the adiabatic curve. 



where lnfc cn d = \J 47r/3r0 cn d/ (mm p ) ~ 1/2 is the fidu- 
cial wavenumber at the end of warm inflation. Note that 
the amplitude of adiabatic perturbations is independent 
of the magnitude of the inflationary potential. In Ap- 
pendix B we show that the amplitude of adiabatic per- 
turbations in warm inflation is independent of the mag- 
nitude of the inflationary potential for all polynomial po- 
tentials. 

Figure 2 shows the amplitude of adiabatic, isocurva- 
ture and tensor perturbations for the n = 2 polynomial. 
In Appendix B we derive results for the more general 
polynomial and exponential cases. We find that the ra- 
tio of the amplitude of perturbations agrees with the in- 
equality expressed by equation ([38]). For the n = 2 case 
the isocurvature and tensor perturbations are strongly 
suppressed below the adiabatic perturbations. This dif- 
ference can be reduced by increasing the magnitude of 
the inflationary potential. 

The spectral indices are related by 

3 3e . , 

n s - 1 = n iso = ~n g = - —, (54) 

and are all red. In Appendix B we derive more results 
for the general case of 4> n potentials. There we show that 
for polynomial potentials the spectral index is always red 
unless n > 4, or for the case of an exponential potential 
(see Appendix B), in which case the spectrum is blue. 
The latter open up the interesting possibility of primor- 
dial black hole production in the early universe p7j. A 



scale invariant spectrum, n s = 1, is found only for the 
ii = 4 model. 

The dissipation ratio, r, and the slow-roll parameter e 
are related by 

'"TS (£)■"" (55 » 

so that the ratio of isocurvature to adiabatic perturba- 
tions, given by e/r 2 , is 3(r7i/r) 2 <C 1. 

Figure 3 shows the spectral indices of the adia- 
batic, isocurvature and tensor spectra as a function of 
wavenumber. In all cases the spectral indices starts as a 
flat spectrum, when e <C r. In the case of n = 4 spectral 
indices are exactly flat while for the exponential potential 
the spectral indices are almost flat. 

VI. NUMERICAL MODELS 

While our analytic results simplify the whole analysis 
of warm inflationary models, it is important to verify the 
range of validity of this analysis and find where the as- 
sumptions used in their derivation are accurate. To test 
this we have numerically evolved the inflationary and ra- 
diation fields through a warm inflationary epoch, and 
calculated the adiabatic and tensor perturbation spec- 
tra. We numerically solve equations (Q) and (Q) for the 
evolution of the field and equation ([|) for the evolution 
of the background universe using the methods of [EOl . 
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The results for the evolution of the energy-density of 
the fields for the case of a polynomial n — 2 potential are 
shown in Figure 1. The magnitude of the inflationary 
and radiation energy-density fields are both normalised 
to unity in units of the potential at the start of inflation. 
We set p-y = p,f, initially, as it is reasonable to assume 
that inflation begins from a radiation-dominated epoch. 
We tested a number of other initial conditions including 
p 1 = Q and found little variation in our final results, as 
the models soon settled down to the same evolution after 
the first few e- folds. 

To parameterise the models we used the more gen- 
eral expression for polynomial inflation from Appendix 
B, and choose parameters so that the amplitude of adia- 
batic perturbations Sh = 2 x f 0~ 5 at the 50 th e-fold, and 
allowed a total of 70 e-folds of warm inflation to elapse. 

Figure 1 verifies that not only the shapes of the curves 
are well approximated, but also the relative magnitudes 
of the radiation and inflationary energy-densities are ac- 
curately reproduced. The main assumptions underlying 
our analytic results are that of slow-roll inflation and 
quasi-stable radiation production (p~ <C T0 2 ). The latter 
assumption breaks down at the start of warm inflation as 
the radiation field initially is redshifted away. But radi- 
ation production by dissipation of the inflaton field soon 
catches up with the redshifting and stability is reached 
after a few e-folds. 

We found accurate agreement between our analytic 
models and the full numerical models for polynomial po- 
tentials in the range 0.5 < n < 4. In most cases the 
agreement was < 1% over most of the range, except near 
the beginning of inflation. In the n = 4 case the en- 
ergy densities tracked each other as expected from our 
discussion in Section IV. We also tested our results for 
the exponential potential discussed in Appendix B, and 
again found similar accuracy, although here again the ra- 
diation field drops faster than the energy density in the 
inflationary field and warm inflation has to end by the 
natural decay of the inflaton, or by some other means. 

We have also computed the amplitude and scale de- 
pendence of the perturbation spectra for adiabatic per- 
turbations directly from equations (|l^) and (|2~l|). These 
are plotted for the case of n = 2 in Figure 2. Again close 
agreement was found, with deviations less than 1% over 
the range of wavenumbers. Again the poorest agreement 
was over the few e-folds at the beginning and end of warm 
inflation. 

Finally we tested the predicted spectral indices for adi- 
abatic perturbations (Fig. 3). Again we found good 
agreement between the analytic and numerical models, 
except for the first and last few e-folds. The sudden 
change in n s — 1 at the beginning of warm inflation 
(N w 70) is real and due to the sharp transition to warm 
inflation when the change in the radiation density, p 7 , is 
important. 



VII. OBSERVATIONAL CONSTRAINTS 

In standard inflation, there are four observable con- 
straints (Sh, A g , n s , n g ) for the three parameters of 
slow- roll inflation, Vq, e, r\, where Vq is the magnitude of 
the inflationary potential. With four constraints there is 
a redundancy in the observations allowing a consistency 
relation j p5"| . This is usually expressed as a relationship 
between the tensor-to-scalar ratio and the slope of the 
tensor spectrum, 

R a = (56) 

This appears to be generally true for slow-roll inflation, 
although recent work on pre-heating at the end of super- 
cooled inflation may effect this result [|s|-[l~lf . 

In warm inflation the dissipation factor, r, becomes an 
active parameter (since it is present at the end of super- 
cooled inflation, but assumed inactive) giving four con- 
straints for four parameters. Hence we do not expect the 
consistency relation of standard inflation to hold in warm 
inflation. This relation can only arise by coincidence. If 
the consistency relation is not observed it may indicate 
a warm inflationary phase. The dissipation factor can be 
estimated from the observables by the relation 

r^-n.^RiAf. ,57) 

To discriminate between warm and standard inflation 
requires a measurement of all four observables. The MAP 
and Planck satellite missions will probe the CMB and 
provide strong constraints on the amplitude and shape 
of the scalar perturbation spectrum. Both MAP and 
Planck have polarisation detectors which are sensitive to 
tensor modes, but Planck has a much broader range of 
frequencies which will aid the subtraction of foreground 
contamination of the CMB signal. In practice the tensor 
modes may be so small as to be unobservable, hence it 
may not be easy to distinguish between the two models. 

Current observations, from the COBE satellite, and 
subsequent airborne detectors (see ]30| ] for a recent anal- 
ysis) , can be used to place constraints on the amplitude of 
scalar perturbations and the slope of their spectra. The 
low COBE multipoles indicate that S H w I0 -10 , giving 
the constraint 

V 1 / 4 » 7.6 x 10- 4 (a 1 /3 er -3 ) i/4 mp) 

m 7.6 x I0 15 (a 1/3 er" 3 ) 1 / 4 GeV. (58) 

If r 3 « ot^ 3 s, these energy scales are pretty much in 
keeping with those of Grand Unified Theories. However 
if r 3> 1 as required for warm inflation we open up the 
interesting possibility that the energy scale required can 
be much lower. 
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VIII. CONCLUSIONS 

In this paper we have investigated the phenomenology 
of the warm inflationary scenario. In the slow-roll, warm 
inflation epoch we have found the general relationship 
between the energy- density in the inflaton and radiation 
fields. This has led us to the general criteria for warm 
inflation, e < r, \r)\ < 3r 2 and r > 1. The usual inflation 
condition, e < I, is generalised as dissipation can help 
sustain accelerated expansion of the universe for a wider 
range of inflationary potentials. 

Generalising the analysis of inflation in the slow-roll 
regime, we have derived expressions for the perturbation 
spectra of adiabatic and isocurvature fluctuations gener- 
ated by thermal fluctuations during the warm inflation- 
ary epoch. The adiabatic spectrum produced by thermal 
fluctuations dominates that produced by quantum fluctu- 
ation. Its amplitude and spectral slope are both different 
from quantum generated perturbations and can lead to 
red or blue spectra. Isocurvature perturbations also are 
generated due to internal fluctuation between radiation 
fields. These are generally suppressed compared to the 
adiabatic spectrum. 

Gravitational waves are still generated via quantum 
excitation, and have the same amplitude but different 
slope from standard inflation. The tensor-to-scalar ra- 
tio is suppressed by the relatively high amplitude of the 
adiabatic spectrum. 

No consistency relation exists in warm inflation as 
there are now four parameters (the standard three of in- 
flation and an active dissipation factor) and in general 
only four observables. This is probably the most im- 
portant observational test of warm inflation, but may be 
difficult, even for Planck, if the tensor modes are strongly 
suppressed. In this case measurement of the isocurvature 
modes relative to adiabatic would provide an alternative 
test. However, verification of the consistency relation 
would rule out substantial interactions during inflation. 

There are two important differences between standard 
and warm inflation. First, the slope of the inflationary 
potential does not need to be small, as dissipation will 
allow slow-roll even for steep potentials. Second, the en- 
ergy scale of warm inflation may be much lower than re- 
quired for standard inflation, due to the higher amplitude 
of thermal adiabatic perturbations and the constraints of 
Sachs- Wolfe effects observed in the CMB. 

To investigate the warm inflationary model we need to 
study further the evolution of scalar perturbations dur- 
ing the transition phase from warm inflation to radiation 
domination when the energy densities of the radiation 
field and inflationary field are comparable. This will oc- 
cur more smoothly than in standard inflation and should 
lead to observational differences between the models. In 
addition the model we have presented here for isocur- 
vature perturbations is crude. A more complete under- 
standing of the generic prediction of isocurvature per- 
turbations for warm inflation is needed before detailed 



predictions can be made. 

The model we have presented here is an extreme one, 
in that we have assumed that it is dissipation-dominated. 
A whole range of models lie between that of warm infla- 
tion and the standard supercooled inflation, and it would 
be interesting to explore their properties in the general 
framework we have outlined here. The development in 
|]f9| also would be beneficial for this. 

Finally these results should be propagated to the post- 
recombination era using general-purpose transfer codes 
such as CMBfast However to make definite, high- 

quality predictions for the CMB, such codes must be gen- 
eralised to remove any of the relations of standard super- 
cooled inflation, such as the consistency relation. This 
will allow us to predict with some confidence whether the 
MAP and Planck missions will be able to probe the high 
redshift universe with sufficient accuracy to distinguish 
between warm and supercooled inflation. 

Acknowledgements: ANT and AB acknowledge the 
support of PPARC, and thank David Wands for useful 
discussion about the evolution of non-adiabatic super- 
horizon perturbations. 

APPENDIX A 

The origin of the "freeze-out" wavenumber ItF in equa- 
tion (|l|) is reviewed here following the original and more 
detailed explanation in Jl4|]. To account for fluctuations, 
the inflaton now is written as <fi(x,t) = <fio(t) + S<fi(x,t). 
Assuming the universe remains near thermal equilibrium 
during warm inflation, the fluctuations of the inflaton 
field obey a fluctuation-dissipation relation. The evo- 
lution equation for the fluctuations, 6(j)(x,t), is obtained 
from the linearized deviation of the zero mode equation of 
motion equation ([!]) , with inclusion of the spatial Lapla- 
cian term and the addition of a white-noise random force 
term 

r d^fc 1 t) = _ [k2 + y„ {(f)o)]Sm t) + e(fe] t)j (59) 

where V" (0o) = d 2 V((j)o)/d(j)Q and we have Fourier trans- 
formed to momentum space. The white-noise represents 
the action of the thermal heat bath on the inflaton field, 
with the spectrum of the this noise term governed by the 
fluctuation-dissipation theorem, 

(f(M)>=0 (60) 

and 

(e(M)£(-fc',0)? T =°° 2TT(2Tr) 3 S^{k-k')S(t-t / ). 

(61) 

The solution of equation (^) is 

S(f>(k,t) » i e -(*-*o>M« / e^WW^fc,^ 

r J to 

+ ^(fc,io)e- (t ~' o)/r(0o) , (62) 
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where r(<p) = Y/(k 2 + V" ((/>)). The first term on the right 
hand side acts to thermalize S<p, whereas the latter is the 
memory term for the initial value of 5(f), which over time 
becomes negligible. 

In the cosmological case, equation ([39]) must be in- 
terpreted for physical wavenumbers, since the thermal 
effects as expressed in this equation act in accordance 
with the physical wavenumbers of the system. Thus for 
a given fluctuation mode, 5<fi(k c ), at comoving wavenum- 
ber k c , equation ( |59"|) expresses its dynamics within the 
time interval 1/H when the physical mode is k p i iys = 
k c e~ Ht w k. To be consistent with the thermalization 
conditions, during the time interval ~ 1/H the mode 
4>(k c ) must thermalize at physical scale k. In terms of the 
solution equation ([32]) , this requires the memory term to 
be negligible within a Hubble time, thus 



fc 2 + v"{<M 

HY 



> 1. 



(63) 



The freeze-out wavenumber kp is at the point where this 
condition first holds, which for V"(<fio) < YH is 

k F = VTH. (64) 

If V"(<po) > HT then the freeze-out wavenumber is the 
same as for supercooled inflation, kp — H. However, in 
general this regime never occurs during warm inflation 
and equation (p4) is the governing condition. Note that 
H varies during warm inflation, thus based on equation 
(|64|), so too will kp. However, the variation of H is typ- 
ically very small, less than a factor 10, during the entire 
warm inflation period. Thus, up to order one factors, 
one can treat kp as a constant with H evaluated at some 
appropriate time during warm inflation, such as the be- 
ginning. In our analytic approach we have included the 
</>-dependence of kp . 

In our treatment here we have taken the dissipative 
coefficient T to be independent of fc, whereas in [|14j it 
was mentioned that in general Y will depend on fc and 
generally will decrease with increasing wavenumber. In 
this case, the condition equation (64) still is valid except 
it must now be solved for kp due to the fc-dependence in 

r. 



APPENDIX B 

Polynomial potentials 

In this appendix we explore the properties of the warm 
inflationary scenario for the family of models with the 
functional form 



V(4>) = Am 4 - 
m 



(65) 



where A is a dimensionless parameter, which allows more 
freedom in the normalisation of the potential. 



The slow-roll parameters are given by 



? 2 

n nip 



167T0 2 ' 

with the relation 



Sue)) 2 



2(n- 1) 



(66) 



(67) 



The expansion of the universe is governed by the Hubble 
parameter 



H{4>) = 



8ttA 



n/2 



(68) 



where the dimensionless dissipation parameter is 



1 



i / 1 p 

rn 



-n/2 



The energy density of the radiation field is given by 

(4-3n)/2 



(69) 



3A 3 n m z 



(70) 



2vr 8 T 

and related to the energy-density of the inflation field by 

2/r, 



3 n m Ti 



" 1 = ^lf^ 



A 



4-n 



2 P0 

mi 



(37.-4) 



(71) 



This last relation tells us that for an n = 4 potential 
the energy densities are proportional, p 7 oc p^, and that 
warm inflation cannot end by the condition given by 
equation (|l4|). When the potential is steeper than this, 
n > 4, or for the exponential potential, the radiation field 
drops more rapidly than the inflation field. Although 
warm inflation occurs, the effects of the radiation field 
will diminish and inflation will begin to evolve in the 
same way as for standard inflation. Hence for polyno- 
mial inflationary potentials we require n < 4 for strong 
warm inflation to take hold. 

The number of e- folds, N(<j>), the universe expands 
during inflation is 



8Y 



3A n(4 — n)m p 



(72) 



where we count the number of e-folds from the start of 
warm inflation (note in Section V we counted e-folds from 
the end of warm inflation) . The number of e-folds can be 
related to e and r by 



2n 



111 

to 



(73) 
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If n < 4 then N is dominated by the ratio of e/r at the 
start of the inflationary phase, and to a good approxima- 
tion the total number of e-folds is 



Nt, 



2n \ tq 



4 — n / Eq 

The adiabatic spectrum of density perturbations is 

1/8 / r g N 1/4 / M x 3(4-n)/8 



(74) 



75 V 3 5 tt 7 A 3 



(75) 



Substituting in equation (|7|) we find that the amplitude 
of adiabatic perturbations in warm inflation is indepen- 
dent of the magnitude of the inflationary potential for 
polynomial potentials, when expressed as a function of 
e-folds, or scale. 

The spectral index of adiabatic fluctuations is related 
to the slow-roll parameter, e, and the dissipation param- 
eter, r, by 



n s — 1 = - 



(4 - n) 3e _ (4 - n) 3 _ 3(4 - n) 
n 47 ~~ n 8™ 9 ~ (5n - 4)" L ' 



(76) 



where the last two expressions relate the adiabatic spec- 
tral index to the isocurvature and tensor spectral index. 
From this expression we see that the warm inflation- 
ary adiabatic spectrum is only exactly scale invariant 
(ri s = 1) for potentials with n = 4, and is red unless 
n > 4. 

The isocurvature spectrum is given by 



6L = 



X'' 



1 



50 V 6 3 tt 15 



1/8 



_r 

a m 



(5ra-4)/2' 



1/4 



(77) 

The ratio of the amplitude of isocurvature to adiabatic 
perturbations is given by 



Rr = 



16irr 2 



32tt 



n-1 



(78) 



In the special case of n = 2 this is a constant, as found 
in Section V. 



r] = 2e = 



8ttk 2 



(80) 



The energy densities of the radiation and inflationary 
fields are related by 



Pi 



_3_ / nip \ 
2vr V8rK 2 / 



3/2 



(81) 



Hence the energy density of the radiation field falls more 
rapidly than that of the inflationary field, as predicted. 
In the exponential model, warm inflation will end in the 
same way as discussed earlier for n > 4 polynomial po- 
tentials. 

The number of e-folds from the start of warm inflation 

is 



3 \ m p 
r\ rr 
e) o \e 



(82) 



The amplitude of adiabatic scalar perturbations is 



8 H = ■==[- 



8 /2\ 7/8 1 



75 V 7T / 3 5 /s a i/4 V m . 



(83) 



while that of isocurvature perturbations is given by 



with ratio 



16%/2 \ m p / \ "'/. 

The adiabatic and isocurvature spectral indices are 

5e 



(85) 



1 3e 
n, — 1 = — . 

4r' 



4r 



(86) 



For warm inflation with an exponential potential, the 
adiabatic spectrum is always blue in contrast with stan- 
dard inflation while the isocurvature is always red. The 
spectral indices are related by 



Exponential potential 



n s — 1 = — -ni a 
5 



(87) 



Consider potentials of the form 

y(0) = ^ o ex P ^, (79) 



K 



where k is a free parameter with units of mass. r , . „ „ „ _ , 
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